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SUMMARY 


Cross sections for the reaction H + + H(ls) -> H(ls) + H + are calc ula ted 
for incident particle energies between 10 ev and 2250 ev using an impact 
parameter method. The required electronic energy levels were found by varying 
one parameter in the molecular ion wave function at every value of proton 
separation. The resulting cross sections are compared with those obtained 
using more accurate wave functions and with experimental data. 


INTRODUCTION 


Interest in charge exchange cross sections has increased in recent years 
because of their usefulness in explaining the presence of certain ions in the 
upper atmosphere and because of their application to ion beam techniques. 

Methods for calculating cross sections were described by Bates, Massey, and 
Stewart (ref. l) . Their "perturbed stationary state" (PSS) method was used by 
Dalgarno and Yadav (ref. 2 ) to calculate cross sections for the resonant reaction 
a + H(ls) -* H(ls) + H at energies too low for use of the Born approximation. 
Exact molecular ion wave functions were used in the calculation. For analogous 
reactions among larger atoms, exact wave functions are not available. Approxi- 
mate wave functions with three variable parameters were used by Moiseiwitsch 
(ref. 3) to calculate cross sections for the reaction He + + He -* He + He + , 
but the parameters were optimized only for infinite separation of the nuclei. 

In order to examine the errors introduced by the use of approximate wave 
functions, calculations of cross sections for the reaction H + + H(ls) -*■ H(ls)+H + 
were performed at the Ames Research Center by means of the impact parameter 
method of reference 1, using hydrogen- like wave functions with a variable param- 
eter. In this method as in the PSS method, the adiabatic approximation is used; 
that is, it is assumed that the nuclei move much more slowly than the electron, 
so that the latter can be described by stationary state wave functions which 
must be determined for each internuclear distance . For a hydrogen molecular ion 
only one parameter is needed to obtain a fairly good approximation to the wave 



function (ref. 4). An algebraic expression for this parameter in terms of the 
particle separation has been obtained by using the variational principle, and 
the resulting energy eigenvalues were used in conjunction with the equations of 
reference 1 to obtain cross sections for the reaction H + + H(ls) -*• H(ls) + H 
for incident particle energies between 10 ev and 2250 ev . The calculations were 
programmed for the IBM 709 computer by Miss Marilyn A. Myers. 


SYMBOLS 


a proton moving in the laboratory frame of reference 

b proton stationary in the laboratory frame of reference 

e base of natural logarithms 

e_ charge on an electron 

E initial kinetic energy of proton a 

E 0 ground state energy of the electron in a hydrogen atom 

■fi Planck's constant divided by 2it 

H Hamiltonian operator 

1 -Fl 

m mass of the electron 

Q charge exchange cross section 

r a distance from electron to proton a 

r , distance from electron to proton b 

b 

R distance between protons a and b 

R c (p) distance of closest approach of protons a and b 

t time 

v(R) relative velocity of protons a and b 

v(°°) relative velocity of protons a and b at infinite separation 

X + symmetric wave function 

X" antisymmetric wave function 

2 



p 


variational parameter in equation (b) 
energy eigenvalue of the symmetric state 
e“ energy eigenvalue of the antisymmetric state 

£(p) phase shift defined in equation ( 3 ) 

p distance of closest approach of protons a and b assuming that their 

paths are undeviated in the encounter 

t wave function of the electron 

ANALYSIS AND DISCUSSION 


Consider a proton a approaching a second proton b at constant speed, 
initially directed along a straight line whose closest distance to proton b is 
p as shown in the sketch. An electron e_ is originally in the Is state of 
system a. Its distance from proton a is r a , from proton b, r^ . Since 


a 



protons a and b are identical, the electronic wave function must be symmetric 
or antisymmetric with respect to interchange of the two protons. Let X + be an 
appropriate symmetric wave function, X” an antisymmetric wave function, and let 
the initial wave function of the system be (l//2) (X + + X )exp[- (i/fi)E 0 t ] . After 
the collision the electronic wave function is (ref. l) : 



X + exp 




E 0 ) dt 


+ X~exp 




(e“ - E 0 )dt 



( 1 ) 


in the approximation that there are no transitions to higher states, where E 0 
is the ground state energy of the electron at infinite separation of the protons, 
e + is the ground state energy of the symmetric state, €" is the ground state 
energy of the antisymmetric state. The cross section for charge exchange is 
(ref. 1) 


3 



Q, = 2jt 


CO 


( 2 ) 


p si 


C(p)<lp 


where 


5(p) 





v 2(o°)p2 

V 2 (R)R 2 _ 


(3) 


where v(r) is the relative velocity of the two protons, and R c (p) is the 
distance of closest approach- Moiseiwitsch (ref. 3) used the further approxi- 
mations v(R) = v(°°) and R c (p) = P which we shall use for the case of hydrogen. 

The symmetric and antisymmetric spatial wave functions used in our 
calculations are, respectively, 


X 





(4) 


where 3 is the variational parameter. Inserting these into the expression for 
the initial wave function shows that the electron was assumed to he originally 
in a Is state near proton a. The energy eigenvalues for the symmetric and 
antisymmetric states are, respectively. 


where 


+ = /X + HX + dr 
/X + X + dT 


and 


- = /X'HX dr 
/X"X" dT 



(5) 


is the Hami ltonian operator for the electron in a hydrogen molecular ion, m is 
the electron mass, e_ is the electron charge, dT is the increment of volume, 
and the integrations are carried out over all space . Evaluation of the integrals 
is simple except where r a and r^ appear together in the integrand. In the 
latter case it is convenient to use confocal elliptic coordinates (ref. 4). The 
integrals of equation (5) are separately minimized with respect to 3 to deter- 
mine e + and e“ in terms of R. The integrations for £(p) were carried out 
on an IBM 709 computer for many values of p. These values were then used in 
the integration for Q at several values of v(°°) . 

A plot of the energy eigenvalues of the symmetric and antisymmetric states 
is shown in figure 1. Plots of the integrand of equation (2) appear in 
figures 2 and 3* The error in assuming Rc(p) = P in equation (3) may he 
estimated by noting that the e + curve determines the classical limit of prox- 
imity of the protons. Contributions to Q for smaller values of p should be 
omitted. However, figures 2 and 3 show that these are very small anyway. In 
fact, the contribution is only 1 percent at 10 ev and gets smaller for higher 
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energies. At higher energies the effect of assuming v(R) constant is mainly 
to ignore the slowing down of the protons when they are very close (p < l). 

Figure 3 shows that for small p the resulting change in £(p) would make little 
difference , since the integrand is very small. This is not the case at low 
energies where the velocities are changing by a significant fraction at several 
Bohr radii. Thus the calculation is not made below 10 ev. Ferguson (ref. 5 ) 
finds that consideration of momentum transfer of the electron (neglected in our 
approximation) reduces the cross section by 5 percent at 2000 ev and by larger 
amounts at higher energies , so our calculation was not extended much above 
2000 ev. 

The resonant charge exchange cross sections obtained using the variational 
wave functions are shown in the following table at various incident particle 
energies E: 


E (ev) 

v(°°)(l0 6 cm per second) 

Q (10 ^cm 2 ) 

10 

4.38 

38.5 

40 

8.76 

30.4 

100 

13.84 

25.6 

160 

17.51 

230 

250 

21.9 

21.3 

1000 

43.8 

15-3 

2250 

65-7 

12.2 


They are compared with those of reference 2 and with the experimental results of 
Stebbings, Smith, and Fite (ref. 6) in figure k . The good agreement between 
these results and those obtained by more exact calculations suggests that the 
variational method may be applied with satisfactory results to more complicated 
atoms . 
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Experimental data 
Ref. 6 






